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I Abstract. Using the multiple scales method, the interaction between two bright and one 



dark solitons is studied. Provided that a long wave-short wave resonance condition is sa- 
tisfied, the two-component Zakharov- Yajima-Oikawa (ZYO) completely integrable system 
is obtained. By using a Madelung fluid description, the one-soliton solutions of the corre- 
sponding ZYO system are determined. Furthermore, a discussion on the interaction between 
one bright and two dark solitons is presented. In particular, this problem is reduced to solve 



^ I a one-component ZYO system in the resonance conditions. 
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1 Introduction 

In many physical applications more than one single wave is propagating in a nonlinear medium, 
^ i and the interaction of several waves has to be taken into account. In particular, to describe the 

. propagation of three nonlinear pulses in some dispersive material one has to solve simultaneously 

a set of coupled nonlinear Schrodinger (NLS) equations. For instance, it occurs in: (i) the 
propagation of solitonlike pulses in birefringent optical fibers [HElElll]; (ii) the nonlinear wave 
dynamics of Bose-Einstein condensates [5l [6] ; (iii) the soliton propagation through optical fiber 
array [U EJ El [9] ; (iv) the nonlinear dynamics of gravity waves in crossing sea states [10] . 

In this paper, a theoretical investigation on the interaction between two bright solitons (short 
waves) and a dark soliton (long wave) is carried out on the basis of previous investigations that 
have reduced, under suitable conditions, the study of bright-dark soliton interaction to the 
study of the existence of a "long wave-short wave resonance" (LW-SW resonance) [11]. This 
resonance phenomena has quite a large universality. For instance, in plasma physics it describes 
Langmuir solitons moving near the speed of sound |12l [T3l I14j . in hydrodynamics it appears 
in the study of internal gravity waves [I^ and a general study of LW-SW resonance [16^ 117] . 



*This paper is a contribution to the Proceedings of the Conference "Symmetries and Integrability of 
Difference Equations (SIDE-9)" (June 14-18, 2010, Varna, Bulgaria). The full collection is available at 
[iittp://www.emis.de/journals/SIGMA/SIDE-9.html| 
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in quasi-one-dimensional molecular crystals it describes the resonance between the exciton and 
phonon fields in Davydov's model \W \ \7U[ I21j. Recent extensions of the LW-SW resonance 
to two dimensions and more components have been discussed and solved by several authors 
[2211231 [21]. 

In the resonance condition, we try to reduce our problem to the study of the completely inte- 
grable system of two-component Zakharov-Yajima-Oikawa [12^ [T3] (see equations (9) and (10) 
of [H]). Besides its relevance in nonlinear optics, the same system describes the Davydov model 
with two excitonic components [25] . It is interesting to mention that a vector-like generalization 
of Zakharov-Benney's equations for long-short wave interaction were obtained 25 years ago in 
the study of magnon-phonon interaction in a many sublattice XY chain [26] , which in the case 
of one-sublattice XY chain reduces to the Yajima-Oikawa system. In the next section, the basic 
equations describing the three wave interaction will be presented and, using a multiple scales 
analysis, the Zakharov-Yajima-Oikawa system is obtained. In Section[3l the Madelung fluid de- 
scription is used to discuss analytically the above system. In particular, several solitary solutions 
in the form of traveling waves are presented. In Section HI the interaction between one-bright 
and two-dark solitons is discussed, and a simple one-component Zakharov-Yajima-Oikawa sys- 
tem is obtained in resonance condition. In particular, one soliton solution is presented. Finally, 
remarks and conclusion are presented in Section [5l 



2 Basic equations and multiple scales analysis 

We consider three nonlinear dispersive waves propagating in an optical fiber. We assume that 
these waves are associated with weakly nonlinear dispersion relations, here denoted by Uj = 
u}j{kj : l^ip, |^3p)) j = 152,3, respectively. Each Aj stands for a complex amplitude 

that, due to the weakly nonlinear dynamics of the medium, is affected by a modulation in both 
space and time. As it is well known [27], due to this dynamics, each mode of the system can 
be represented by a wave packet where a carrier wave is amplitude modulated. In order to 
obtain the evolution equation for each mode, we apply the well-known method of Karpman and 
Kruskal |28| (see also [271129]). To this end, let us denote by e*('=o^-'^o*) the basic carrier wave. 
Then a Taylor expansion around {ko,uJo) and = of each will give 
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Replacing — — —i-§t, h — ko — i-^, after a translation of coordinate (x — )■ x — (f||)o*); 
the following nonlinear system of three interacting waves is obtained 

■ dAi .^^ dAi ai d'^Ai , , ,9 , 1 ^ 19 ^ 1 ^ 19 ^ 

+ ^^1^ + -^^^ + a2\AiYAi + aalAalMi + aA\A^YAi = 0, 



dt dx 2 dx"^ 

dt 2 5x2 



i^r + ^^+ 72|^lP^3 + 73|^2P^3 + 74|^3|'^3 = 0. (1) 



Here we denoted Vi = {^-^ — (^f) > ^ = 1,2 and the constants ai, /3i, 71 are related to 
derivatives of cjj with respect to ki (ex. ai = — i^^r] , ■) while 02, ... ,74 to the derivatives 



with respect to \Ai\'^ (ex. 02 = (^ gfXi]^ ^ ' ' ')■ 
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Further on we shall consider channel 3 with normal dispersion and 1 and 2 with the anomalous 
one [T]. Then following Kivshar [TT] it is convenient to introduce new field variables 

(uq, a{x,t), (j){x,t) being real quantities). Note that for solitary wave solutions the additional 
conditions of vanishing of a at the infinity has to be imposed. Then the equations of Ai and A2 
become, respectively 

i-^ + ^^i-g^ + Y~dx^ ^ ("2|*ir + 031^2! ) ^1 + 2a4Uoa^i + a^a^^i = 0, 

+ ^^2^ + + (/52|1'iP + /33|^2n ^2 + 2/34noa*2 + /34a''f 2 = 0. (2) 

As concerns the ^3-equation, separating the real and the imaginary part, the following system 
of coupled equations is obtained 

da 71 d^4> , , , 

— I uoTT-^ + (nonlmear terms) = 0, 

ot 2 ox"^ 

dcj) 7i d'^a / I ,9 I i2\ / 1 

— — + 274^00 + '^'q^ + (72IV1I + ^31^21 ) + (nonlinear terms) = 0. 

In both these equations the parenthesis (• • • ) contains all the other nonlinear terms which will 
be irrelevant in a multiple scales analysis. From these last two equations the following equation 
satisfied by a{x, t) is easily obtained 

d'^a 2^2^ 2 ^4^ 71 92 . u., |2^ |,T, |2\ /ON 

5^+7i74.o^ + Y5^+-oy^(72|^i| +73|^2|) (3) 
+ (higher order nonlinear terms in (a, (p) and their derivatives) = 0. 

The linear part of the a equation corresponds to an acoustic field with dispersion relation 
(71 < 0, 74 > 0) 



and phase velocity c = Lo/k, where = |7i|74tio- 

We shall perform a multiple scales analysis of the system ([2]) + (l3|) [II] . We introduce new 
scaled variables 

t ^ et, X ^ — ct) 

and new functions 

3 3 

a^ea, (p ^ ecp, ^1 =^ e4^'i, ^2 ^ e*^2- 
Then in order | in e from a equation we obtain 

-2c-^ + uoYo^ (72|^i|' + 73|^2|') = 0. (4) 

All the other terms in a equation contribute to higher order in e. In the order | from ^'j 
equations we obtain Vi = V2 = c. This is the well known long wave-short wave (LW-SW) 
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resonance condition: "t/ie group velocity V of the SW is equal to the phase velocity of the 
LW^ |16j . In the next order (|) in e from the 'I' equations we get 

Ml ai 52^-1 ^ 

^^ + Y^ + ^^^^°"^^ = °- 

The equations ([1]) + (l5|) represent an 1-dimensional 2-components Zakharov [12], Yajima- 
Oikawa [13] system. As mentioned in the Introduction the same system in the same LW-SW 
resonance condition was obtained in a Davydov model with two excitonic modes coupled with 
a phonon field [25]. The same line of reasoning was used in [23] for three interacting waves in 
2-dimensions. 



3 Madelung fluid description 

The special case (oj = ft, 72 = 73) is completely integrable [23] and will be considered in the 
following. In this case, simplifying the notations, the system (j3]) + ([5]) is written in the following 
form (7 > 0, /3 > 0) 

^-7^(i^ip+i^.n=o, 

The equations will be transformed using a Madelung fluid description [30l [31] . We write 

where pi, 9i are real functions of {x,t) and moreover pi are positive quantities. Introducing this 
expression into a-equation this becomes 

da d , . 

__,_(„+p,).0, (6) 

while from the equations, after the separation of real and imaginary parts, we obtain 

dpi d , , 
^ + ^(-.A) = 0, 

which is a continuity equation for the fluid densities pi = \^%\' with Vi{x,t) = ^^g^'*^ the fluid 
velocities components and 

_^ + l_L^_ir^V-,3a = 0. (7) 



dt 2 ^/pi dx"^ 2 \ dx J 

Differentiating this last expression with respect to x the following equations of motion for the 
fluid velocities Vi are obtained 

9 ^ ^ 9 \ _ 1 g / 1 ^da 
dt ^ dx J ^ 2dx \ y/pi dx'^ J dx ' 

In the right hand side of ([8]) a{x, t) plays the role of external potential, and the first term is the 
derivative of the so called Bohm potential, 1^%?^, and contains all the diffraction effects 
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(quantum effects in quantum problems). By a series of transformations the equation ([8]) is 
written as |31| 



dvi dpi 



^ + 1^ _^ .Aq _ 2/3a— = (9) 
dx 4 dx^ ^ dx 9x ' 



where Cj are arbitrary integration constants with respect to x, eventually time dependent. Al- 
though ([9]) seems to be more complicated then the initial equation ([7]), it can be solved in two 
special situations, namely 

• motion with constant velocities vi = V2 = vq, 

• motion with stationary profile current velocity, when all the quantities pi{x,t), Vi{x,t), 
a{x,t) are depending on x and t through the combination ^ = x — uqI. 

Both cases will be analyzed in the following. 

3.1 Motion with constant velocity (v^ = ^2 = Vq) 

In this case from the continuity equations ([6]) one sees that both pi{x,t) and p2{x,t) depend on 
^ = X — vot. We assume that also a(x, t) depends only on ^. Then the a equation gives 

a = -p{pi + P2), = — (10) 

and the equations ([9]) write 

where by Ei we denoted — (2cj — Vq), 13 p — )■ p. We shall discuss firstly the situation Ei = E2, 
the discussion of the more general case Ei ^ E2 being postponed for the next subsection. Then 
the equations PT|) becomes 

These are exactly the equations obtained in the case of Manakov's model [32] and extensively 
discussed by us in [331 [MJ. In the following we shall present several periodic and traveling wave 
solutions of (fT2]) . 

It is convenient to introduce the quantities z+ = pi + p2) and = pi — p2', they satisfy the 
following equations 

ld^Z4- „dz+ 3 dz"^ 



- + -p^ = 0, 



4 d^ 2^ dC 

Id^z^ ^dz^ dz4. dz- , , 

The second equation is a linear differential equation for z_ once z+ is known. A special solution 
is 

= (Pl - P2)^+> Pi + P2 = 1> 

which together with the definition of 2+ gives 
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and the problem is reduced in this simple case to find a solution of z_|_-equation. This integrated 
twice gives 



1 / dz4 ^ ^ 



-^lzl + Ezl + Az+ + B = P3{z+), (14) 



where ^3(2+) is a third order polynomial in z^. The periodic solution of this equation are easily 
expressed through Jacobi elliptic functions. 

For constant velocities, as it is easily seen from ([5]), the densities pi have to satisfy the 
additional conditions 

1 1 d'^ 



2 dx"^ 

which for the previous solutions are satisfied X = E/2 and B = 0. 

Now let us assume that the third order polynomial P3(z_|_) has three distinct roots 

P3(z+) = -l^{z+ - zi){z+ - Z2){z+ - za). 

The restriction B = Q means that one of the roots Z2 or Z3 is zero. We are interested in positive 
solutions of (fT3|) for which ^3(2+) is also positive. The periodic solutions of can be expressed 
through Jacobi elliptic functions and taking into account the positivity requirement mentioned 
before we identify two acceptable situations [35] 

zi > 0, Z2 = 0, Z3 < 0, z+ = zi cr?u, (15) 



9 Zi + l^sl ^zi + |Z3| 

and 

Z3 = 0, < Z2 < ^1, Zj^ = Zi — {zi — Z2)si^U, (16) 

n = ^C, k'='-^, g- 2 



9 zi ^zi 

Solitary wave solutions are obtained in the limiting case k = 1 when cnn — t- sechn, snn 
tanh u, and both solutions (fT5]) and (fT6|) become a bright soliton 

1 2^ 2 

z+^zi —, 



cosh u g ^Jz\ 

It is clear that in this case no energy transfer between the two components takes place. 
The phase 6{x,t) is easily calculated writing 9i{x,t) = vqx + 7j(t); then using ^ we get 

9i = Vox - Qwg - * + ^i- 
As far as the field a is concerned, the solitary wave solution is 

a{u) = -2— 

cosh u 

and, for pzi < uq, the field ^3 describes a grey solution (the inequality gives the dark solution). 
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3.2 Motion with stationary-profile current velocity 

In the case when all the functions depend only on ^ = x — uot integrating the continuity equa- 
tion ([6]) we get 

Vi{x,t) = no + — , 

Pi 

with Ai some integration constants. It is easily seen that the equations of motion keep the same 
form as ([U]) with Ei = -{2ci + Uq). 

For El = E2, the results are almost the same as those of the previous subsection. 

However, for Ei ^ E2, the system (llip will be solved by using a direct method. To this end, 
we look for solutions of the form 

Pi = Ai + Bi sn n, u = 2A^, 

with Ai, Bi, A constants to be determined. Introducing into (jlip we get 



i?i + i?2 = -^, (17) 

-[4A2(1 + + Ei]Bi + + B2)Ai + 2^(Ai + as)^^ = 0. 
Defining the following new quantities Cj, 6j, cq and 5, i.e., 

4A2^ _ AX^k^ 

Bi — Oj, Ai — fflj, 

^1 = 4A2fc2(eo + 6), E2 = 4X^k^{eo -6), S>0 
the first equation (fT7|) gives 
&i + &2 = 1, 

while from the second, after a little algebra, we get 
oi = 3 ( eo + —p— + -5 + 45(1 -61) j 61, 

a2 = ^(eo + ^^-5-45(l-62))62. (18) 

As it is expected this result verify the symmetry condition 1 -H- 2 if J -H- —5. 

Several restrictions result from the positiveness of pi . If both bi are positive quantities smaller 
than unity this requirement implies 

ai>bi> 0. (19) 

Introducing the notation 

the condition (jl9p is satisfied if p > 1. In the limiting case k'^ = 1 the solutions are 

, , 2 ^ 4A2 , , 2 ^ 

pi = (ai — oitanh n), P2 = (02 — 02tanh u) 

p p 

representing shifted bright solitons. 

We can now calculate the field a from (jlOp and using the limit a{u) — )■ as u — )• 00, we get 
ai + 02 = 1. Furthermore, making use of (llSp . this condition becomes 

1 f 1 
where we have used the condition 5i + 62 = 1- 
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4 One bright-two dark solitons interaction 

The previous discussion is easily extended to the situation when two waves have normal dis- 
persion and one is anomalous, corresponding to one bright-two dark solitons interaction. Such 
vector solitons of mixed bright-dark types are of interest in quasi-one-dimensional Bose-Einstein 
condensates (see [6] and references therein). We are interested in a SW-LW resonance regime, 
so several restrictions on the model of three interacting waves will be imposed. With the nota- 
tions of section 2 we shall assume ^ J — ^ Sfca J smd after the corresponding translation of 
coordinates instead of system ([T]) we get 

dAi dAi a\ d'^A\ , . ,o . , . ,9 . , , ,9 . 
'~dr ^ ''^^ ^ Y~d^ ^ a2\Ai\^Ai + a3|^2r^i + a^lAsl^Ai = 0, 

+ f ^ + /32|^i|^^2 + m2\'A2 + PMsl'As = 0, 

i-Q^ + Y~d^ + 72|^lP^3 + 73|^2|% + 74|^3|% = 0, 

where v = (^^^^ — (^^f) ^^^'^ the rest of the notations are the same as in Section [2j 

Further on we shall consider the channels 2 and 3 with normal dispersion and the channel 1 
with anomalous one and introduce new field variables by 

Ai = ^e'^\ A2 = {U2 + a2(x, t))e^(r2*+<A2(x,t))^ ^ ^ ^^(^^ i)y{W4>,{x,t)) 

with U2-, U3, a2{x,t), a3{x,t), (p2ix,t), 4)z{x,t) real quantities, and 

5 = a-iul + a^ul, T2 = l^zul + /34U3, Fs = 73U2 + -fiul- 

The -equation transforms into 

9^- 01^2* , ,2 , X / 2 2^ 

^ ^'"'dx ^ T 9x2" + "2|^r^ + 2(a3'U2a2 + a^uza^i)^ + (0302 + ona^)^ = 0. (20) 

To get the relevant equations for 02, 03 we have proceed like in Section [21 obtaining the following 
equations for the amplitudes 02, 03, i.e., 

-W + y"^^^ + ^^^^ r + ^^^^ J + T + = °' 

9^03 71 52|^|2 / a2a2 d'^a^\ 7? 9^ 03 , , 

^ + T"^^^^ + ^^^^ 9^ + ^^"^9^ J + T 9^ + = °' (''^ 

where, as mentioned before the parenthesis (• • • ) of each equation group all the nonlinear terms 
which will be irrelevant in a multiple scales analysis. Furthermore we shall consider only the 
special case 1*2 = '"3 = ^^O) ^2 = a3, /32 = /33, 72 = 73 and neglect the last nonlinear terms in (pO]) 
and all the parenthesis in (121 p . remaining with the following system of coupled equations 

^~dt ^ ^'"'dx ^ T ar2" + "2|^r^ + 2noa3(a2 + 03)^' = 0, 

-W + y^°^^ + /3i/33no^(a2 + «3) + = 0, (22) 

9202^71 2d\ ^ .^lld'a^ . 



Solutions of Zakharov-Yajima-Oikawa System 



9 



As in ^'-equation appears the combination = a2 + as it is convenient to add the 02- and 
as-equations. In a multiple scales analysis the fourth order derivatives will contribute to higher 
order so we can drop them at the present stage. Then the system (\22\i reduces to 

+ + WITT + a2^r^ + 2uoa3a+^ = 0, 

at ox 2 ox-^ 

-Q^ + y(/5i +71)^ +ng(/3i/33 +7173)^ = 0. (23) 
In the resonance condition (/5i < 0, 71 < 0) 

v = c, = no(|/3i|/33 + I71I73) 

and by using the multiple scales analysis 

3 

x^^/e{x-ct), t^et, a+ — ;> ea+, ^ ^ e^^ 
the system ()23p becomes the well known one-component ZYO system 

+ T 9^ + 2noa3a+^' = 0, 

The one-soliton solution of (j24p is given by 

a+ = Msech^e, * = iVe^'^sech 6, 9 = n{x - xq) - i^t, (j) = ^x - rjt, 

where 

The integration of the equation for a_ is straightforward, and it gives 



M-sech^ei, M_ = -. 



^ (|/3i| - |7i|) iV + (|/3i|/33 - I71I73) m' 



5 Conclusions and remarks 

In this paper, an investigation on the interaction between two bright solitons (anomalous dis- 
persion) and one dark soliton (normal dispersion) has been carried out. The problem has been 
reduced to study a two-component one-dimensional Zakharov-Yajima-Oikawa system in the 
case of long wave-short wave resonance, by using the multiple scales analysis. The system cor- 
responding to the bright solitons has been discussed further using a Madelung fluid description. 
In particular, periodic solutions expressed through Jacobi elliptic functions and stationary so- 
lutions obtained when k'^ = 1 have been presented in two simplifying conditions, namely for 
constant velocity and for motion with stationary profile. 

Remarkably, the above two-component Zakharov-Yajima-Oikawa system is completely in- 
tegrable and multi-solitons solutions can be found using different methods, as the bi-linear 
method of Hirota [23]. However, the Madelung fluid description seems to be useful to find 
various solutions of a generalized ZYO system, containing additional nonlinear terms. 

Additionally, the case of two dark and one bright solitons has been also discussed. Here, 
the use of the multiple scales analysis has allowed us to reduced our problem to a simple one- 
component Zakharov-Yajima-Oikawa system in the resonance condition. In particular, the 
one-soliton solution has been found. 
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Finally, we would like to point out that the discussion for arbitrary values of the parameters 
is still open and investigations in this direction are in progress. Also numerical solutions of the 
dynamical model of three interacting waves, discussed in the present paper, could reveal several 
new aspects of the nonlinearity influence on the system behavior. 
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